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Optical photometric and spectroscopic data are presented that show an asso- 
ciation of an ultraluminous X-ray source in NGC 7331 with a young star cluster 
of mass M = (1.1 ± 0.2) x 10^ Mq and age = 4.25 ± 0.25 Myr. If the ULX is 
part of the bright stellar cluster, then the mass of the progenitor of the compact 
accretor must have been >40-50 Mq in order to already have evolved through 
the supernova stage to a compact object. The companion star is also likely an 
^ ' evolved massive star. The emission line spectrum of the nebula surrounding the 

^ ' cluster can be interpreted as a result of photoionization by the cluster OB stars 

with an additional source of shock excitation producing strong [S ll], [O l] and 
^ N 11 lines. This additional source appears to be as much as five times more 

' powerful than the supernovae and stellar winds in the cluster can provide. Addi- 

^ ■ tional mechanical energy input associated with the ULX itself can help explain 

CSl ■ the residual shock excited line luminosities of the emission region. 

I> 

O I Subject headings: galaxies: individual (NGC 7331) — galaxies: star clusters — 

Q ! line: formation — X-rays: galaxies — X-rays: binaries — X-rays: individual 



X 



(CXOU J223706) 



1. Introduction 

Ultraluminous X-ray sources (ULXs) are point-like non-nuclear sources in external 
galaxies with apparent X-ray luminosities greater than the Eddington limit of stellar mass 
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The existence of ULXs has been known since the first Ein- 



black holes, i.e., Lx^lO^^ ergs s~^ 

stein observations of nearby galaxies (ILong fc van Speybroecklll983l : lHelfandlll984j : iFabbiano 
19891 ) but only with the sub-arcsecond imaging capability of the Chandra X-ray Observa- 
tory have investigations at other wavelengths of the crowded fields of nearby galaxies become 
viable. 

Multi-wavelength investigations enable a much broader range of diagnostic techniques 
in the study of ULXs than can be achieved by X-ray data analysis alone. In particular, 
traditional methods in optical imaging and spectroscopy can be used to constrain the age, 
metallicity and physical state of the local environments of ULXs. This, in turn, can help to 
determine the origin and history of the ULXs themselves and lead to a better understanding 
of the influence of ULXs on their surroundings. 

Many of t he models for ULX s favor young high-mass X-ray binary systems (see, e.g., 
the reviews by IFabbiano fc White! (120061 ): iFabbiand (120061 )). These objects should be pref- 
erentially associated with massive star forming regions such as OB associations and, at the 
high-mass limit, superstar clusters. These environments should be strong sources of line emis- 
sion through photoionization by massive stars, supernovae (SNe) and by the ULX itself and 
may display enhanced forbidden line emission from shock excitation by stellar winds, SNe, 
and by energetic outflows and jets that could acc ompany the ULX X-ray emissi on. Nebulae 



such as those assoc iated with the ULX s in Ho IX (iGrise et al.ll2006a[). NGC 5408 (jSoria et al. 



2005|), NGC 5204 (IGoad et al.ll2002h . and Ho II (IKaaret et al 



2004; Lehmann et al 



200 



J 



are among the best-studied examples of this intricate ULX/star-formation connection. 

As part of a program of spectroscopic studies of the local envrironments of ULX can- 
didates, we have obtained a moderate-resolution long-slit spectrum of a ULX in NGC 7331. 
The X-ray characteristics of this ULX and the optical photometric properties of its environ- 
ment are describ ed in S [21 T he ULX lies roughly in the center of an extended H ll region, 
cataloged as P98 (|Petitlll998l ). and within 30 pc of the center of a bright blue emission region. 
The optical spectrum of this emission region is presented in § [31 We model its spectrum as 
emission from a young star cluster combined with a nebula photoionized by a cluster hosting 
an X-ray source (§ [1]). This model reproduces the observed 4000 — 7300 A continuum spec- 
trum but under-predicts several strong observed emission lines. We argue (§ 14. 3p these lines 
are due, at least in part, to shock excitation associated with the ULX and that radiative 
and mechanical feedback plays an important part in the overall energetics of the emission 
region. We briefly summarize our results and compare the properties of the environment of 
the ULX in NGC 7331 to other ULX nebulae in § E 
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2. X-ray and Optical Observations 



NGC 7331 is a highly-inclined, i = 77° (iGarcia-Gomez et al.ll2002r). SAfslb spiral ga laxy 
subtending 10/5 on the sky. We adopt a distance of 14.6 Mpc (IFerrarese et al.l l2000l ) : 1" 
corresponds to 71 pc. 

Swartz et al. (2004) identified three ULX candidates within the D25 isophote of the 
Chandra field of NGC 7331. All three lie in the northeast portion of the spiral arm struc- 
ture beyond the l arge ring seen in the Spitzer Infrared Nearby Galaxies Survey images 
(IRegan et al.l |2005| ) . An Ha image of this region is shown in Figure [1] with the approximate 
positions of the three ULX candidates superposed. Only the brightest ULX candidate has an 
obvious association with optical emission features in this or any other images we have stud- 
ied. For the remainder of this work, we will refer to this ULX, which is formally designated 
CXOU J223706.6+342620, as CXOU J223706. 



2.1. Chandra observations 

Chandra obtained a 29.5 ks observation of a portion of NGC 7331 on 12 Aug 2001 (ObslD 
2198). The observation was taken in sub-array mode using the Advanced CCD Imaging Spec- 
trometer (ACIS) with the galaxy center imaged near the aimpoint on the back-illuminated 
S3 CCD. CXOU J223706 is approximately 1.'85 off-axis in this observation. Using the known 
position of the nucleu^of NGC 7331, we estimate the absolute positional uncertainty of the 
ULX is ~ 0'.'4. Within a 3o" = 1'.'86 circle about the position of the source, 175 X-ray events 
were detected corresponding to a S/N of 10.6 (in the 0.5-8.0 keV band). Acceptable fits to 
the spectrum can be made using a variety of models. The best-fitting model attempted is 
an absorbed power law (Nh = 6.9l|^ x 10^^ cm'^, F = 2.26l|]:|, = 5.1 for 5 dof after 
grouping to ensure at least 20 counts per spectral bin). The corresponding in trinsic source 



lumin osity is Lx = (2.8±0.6) x 10^^ ergs s ^ according to our previous analysis (jSwartz et al. 



20041 ). (Here and elsewhere, errors are extremes on the single interesting parameter la con- 
fidence intervals.) The luminosity of CXOU J223706 lies at the median value of the 97 ULX 
candidates in our sample of ULXs in spiral galaxies. The average luminosity for that sub- 
sample is 5.9 X 10^^ ergs s~^ and the highest luminosities exceeded several lO'*" ergs s~^. If 
we assume the X-ray luminosity approximates the bolometric luminosity and that the source 
radiates isotropically at the Eddington limit during the Chandra observation, then the mass 
of the compact accretor is 22±5 Mq. 



22/137^04^.096, S = +34°24'56".29 (.T2000): [Arevle fc Clemeirt] (|l99nh . 
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2.2. Hubble observations 



Hubble/WFPC2 observation^ of NGC7331 were carried out on 1997 August 13 and 
14 using F450W an d F814W filters. Astrometric corrections were made using USNO stars 



( iMonet et al.ll2003l ). Tlie la nominal absolute positional uncertainities in the Hubble im- 
ages are (f!3 giving a combined uncertainty of about O'.'S between the Hubble and Chandra 
registration. 

Figure El displays an 8"x8" portion of the F450W Hubble image around CXOU J223706. 
The Hubble image is dominated by two bright extended sources; presumably young stellar 
clusters. The X-ray source is located within the NW cluster, within the formal error. The 
optical source has a FWHM of ~ O'.'l or 7 pc at the distance of NGC 7331 (the source 
is not completely symmetric, Hubble images reveal additional emission in the wing). The 
total background-corrected brightness within a I'.'S diameter circle (equivalent to the seeing 
conditions in the ground-based data, see below) centered on the X-ray source position is 
'^F45ow = 20.44 mag. This corresponds to an absolute magnitud e M45nw = — 10.4 mag. 



lying at about the maximal luminosity of the Galactic open clusters ( lLarsenll200ll ). However, 
taking into account interstellar absorption (§ |4]), the object intrinsic brightness may be 
significantly higher. 



2.3. Isaac Newton Telescope observations 

Archival optical images of NGC 7331, obtained with the Wide Field Camera (WFC) 
on the Isaac Newton Telescope (INT) at La Palmfl, were examined for structure local to 
CXOU J223706. We used images taken in narrow filters [O ill] (central wavelength 5008 A, 
FWHM=100 A), Ha (6568 A, 95 A) and in the Harris broad-band filters B (4298 A, 1065 A), 
V (5425 A, 975 A) and R (6380 A, 1520 A). All the imaging observations were carried out in 
2000 and 2001. We have estimated seeing from the FWHM of unsaturated point-like sources 
in the images to be 0'.'7 in Ha, 0'.'8 in [O ill] and from 0'.'9 to 1'.'2 in the broad bands. 



Astrometric corrections were made using USNO-Bl and USN0-B2 standards ( IMonet et al 



20031 ). The astrometrical error obtained in all these images (using 22 to 27 astrometrical 
standard stars) is less than 0'.'2. We made photometric calibration of the B- and V-band 
images using 19 USNO standard stars that range in brightness from 13.0 to 20.5 mag. A 
systematic error of the USNO photometric standards as high as 0.2-0.3 magnitude is possible 



^Available from ; http: / / archive.stsci.edu/] 
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for the fainter stars (IMonet et al.l 120031 ) . 

We created new Ha and [O ill] images free of continuum emission by dividing the narrow 
band images over the R and V images, respectively. In the resulting images the majority 
of single stars (presumably non-emission-line stars) disappear. From these images one may 
estimate the spatial structure in the emission lines and measure the corresponding emission 
line equivalent widths. Isophotes at the 20, 50 and 85% levels above the background in these 
Ho; and [O III] images are shown in Figure [2] (100% corresponds to the maximum value in 
the 8"x8" field). Taking the 50% isophote as representative of the extent of the emission 
in these lines, one obtains the size of the emission-line region to be ~2" in Ha and [O ill] 
corresponding to 140 pc (FWHM). The peak of the Ha emission is coincident with the NW 
cluster resolved in the Hubble image. The peak [O ill] emission lies between the two clusters 
suggesting each cluster contributes equally to the [O ill] emission. 

The estimated B-band magnitude within a circle of 1'.'5 diameter around the NW cluster 
is IQ'^QTiO.Oe (USN0-B2) and 20'?32±0.06 (USNO-Bl). For modelling the optical spectrum 
(see below), we adopt the value B=20™0. 



3. Optical Spectroscopy 

Optical s pectral data were obtained u sing the 6m Russian telescope with the SCORPIO 



focal reducer (lAfanasiev fc MoiseevI l2005l ) in the long-slit (LS) mode with a O'.'75-wide slit, 
providing spectral resolution of about 10 A. The data cover the spectral range 4000 — 7300 
A. The object was observed on August 29, 2005 at 1'.'4 seeing conditions. The reduction 
process includes all the standard procedures. 

The spectrum is shown in Figure |3l Many strong emission lines are present above a 
well-defined blue continuum. Table 1 lists the observed and dereddened integral fluxes in all 
the emission lines detected. 

We extracted sht images in the strongest lines (H/?, [O ill]. Ha, and [S ll]). The extent 
of the [O III] and Ha emission in the slit images is consistent with the structure seen in 
the INT images. The emission in the slit images can best be described as the sum of a 
point-like structure (FWHM~ 1'.'5 consistent with the seeing conditions) and more extended 
(FWHM~ 5'.'0) emission. That is, the two clusters are not resolved. The centroids of the 
point-like emission in all the lines are coincident within their errors; but in the extended 
component the [O III] centroid is shifted about 0'.'4 to the SSE relative to Ha along the slit 
(cf. Figure [2]). The [S ll] centroid is shifted in the opposite direction, about 0'.'2 to the NNW 
relative to the Ha peak. 
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4. Spectrum modelling 



The size and luminosity of the bright sources in the Hubble blue (F450W) image suggest 
they are compact young star clusters. The INT narrow-band images and SAO slit images 
suggest the star clusters are embedded in a larger emission line nebula (with the slit sampling 
most of the cluster emission but not all of the larger emission line region). We constructed 
synthetic spectra based on a model for the emission consisting of an underlying young stellar 
population and a photoionized nebula. Starburst99 (IVazquez fc Leithereiil2005l ) was used for 
the stellar contribution which provides model spectra for the continuum emission as functions 
of age and metallici ty of the cluster and normalized by the cluster mass. CLOUDY96.01 
( iFerland et al.lll998[ ) was used to simulate the emission line spectrum from the photoionized 
nebula as a function of the input photoionizing source spectrum and ionization parameter, 
the nebular abundances, and the density. The spectral modelling began by finding the best- 
fitting cluster model. This model was then subtracted from the observed spectrum and the 
best-fitting model for the residual spectrum was found using the photoionization code. 



4.1. The contribution from the Star Cluster 



The oxygen abundance and abundance gradient in NGC 7331 is given by lPilyugin et al. 



( 120031 ). At the location of CXOU J223706, about 90" from the center of the galaxy, 12 + 
log(0/H) = 8.32 or about 2-3 times less than the solar value. Using this as a guide, 
we constructed Starburst99 models with [Fe/H]=0.0, —0.4, —0.7. We computed a grid 
of models for cluster ages, tc, ranging from 1 to 30 Myr. The model spectra were fit to 
the observed spect ra allowing for a var iable interstellar absorption component (using the 
reddening curve of ICardelli et al.l (Il998l )). Since the Starburst99 stellar cluster model does 
not predict a nebular emission line spectrum it was necessary to exclude all the bright 
observed emission lines in the model fitting. We find the best-fit parameters of the cluster 
model are tc = (4.25 ± 0.25) Myr, Ay = l'°43 ± 0'P05, M = (1.1 ± 0.2) x 10^ M© where 
M is the cluster mass which scales linearly with the overall model normalization. The Xr 
statistic for this model is 1.30 for ~200 degrees of freedom (the spectral range divided by 
the spectral resolution element with emission lines omitted and reduced by the number of 
model fit parameters). Solar metallicity provides a significantly better fit. The best-fit model 
spectrum (scaled by a factor of 1/3 for clarity) is shown as the lower curve in Figure O 

The age and metallicity derived here are consistent with the presence of the broad 
but weak emission features observed at 4650 A and 5806 A. We interpret these features 
as th e C lll-|-N in and C IV line blends, respe ctively, due to the presence of Wolf-Rayet 
stars (IConti et al.lll983l : ISchaerer fc Vaccalll998l ) in the cluster. We estimate the equivalent 
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widths to be 1.7±0.2 A for the blue blend and 4.2±0.2 A for the red blend of lines. (The 
quoted errors are formal errors using a piece-wise linear approximation to the local continua; 
inspection of Figure [5] suggests the uncertainties are likely higher.) The theoretical blue-to- 
red ratio (ISchaerer fc Vaccalll998l ) most closely approaches the observed ratio at 3-6 Myr 
where the theoretical ratio is ~1 for a wide range of metallicities. 



4.2. The contribution from photoionization 

The best-fit stellar population model spectrum was subtracted from the dereddened 
observed spectrum to give an estimate of the purely nebular spectrum. We used the 3 < 
tc < Q Myr UV spectrum from the best-fit unabsorbed cluster model plus an additional X-ray 
component based on the Chandra spectrum of CXOU J223706 (F = 2.2, 3.2x10^^ photons 

blueward of the Ly edge) for the photoionization source for the nebula model. We 
calculated a set of spherically-symmetric photoionization models with fixed inner (1 pc) and 
outer (30 pc) radii using CLOUDY 96.01. The inner radius is typical for a massive young star 



cluster (e.g., iNilakshi et al.l (120021 )). The outer radius corresponds to the half- width of the 
slit in projection (at the galaxy distance of 14.6 Mpc). We have found that the choice of inner 
and outer radii do not strongly affect the results. The nebula gas density and the low-energy 
cutoff for the ionizing X-ray spectrum were allowed to vary in the modelling. The best- 
fitting gas density is about 30 cm~^ (consistent with the observed [S ll] A6717/A6731 ratio). 
The models proved insensitive to the additional X-ray and corresponding EUV radiation 
from the ULX. This was anticipated because the number flux from the cluster model in the 
Ly continuum (~10^^ s"^) far exceeds the flux of X-ray photons even after extrapolating this 
rather steep X-ray spectrum to the Ly edge. 

The results of the photoionization model predictions and the observed emission line 
parameters are listed in Table 1. For each emission line observed, the ratio of the observed 
(Column 2) and dereddened (Column 3; Ay = 1™4) line fluxes to the H/5 flux are tabulated 
along with the integral dereddened line luminosity (Column 4), the photoionization model 
predicted line luminosity (Column 5), and the residual line luminosity after subtracting 
the photoionization model prediction (Column 6). The best-fit stellar continuum model 
combined with the best-fit photoionization model (convolved with the instrumental response 
function) is shown as the middle (dotted) curve in Figure [3] (displaced by 1/2 for clarity). 
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4.3. The contribution from shocks 



These pure photoionization models predict strong [O ill], He I, and H Balmer Une 
emission excited by the UV emission from the underlying star cluster but are unable to 
reproduce the observed strength of several key emission lines (see Table 1 and Figure [3]). 
A significant excess is observed in [S ll] AA6716, 6731, [Ol] AA6300, 6364 and [Nil] AA6548, 
6583 A doublets. We find the additional dereddened luminosity in the [S ll] doublet, for 
example, is L[sii] = 1.72 x 10^^ ergs s'^. Further, the residual [O ill] A(4959 + 5007)/A4363 
ratio is low, ~20, compared to typical photoionized H ll regions. 

The observed excesses in these low-excita tion lines and the l ow [O ill] ratio are likely 
caused by shock excitation in the nebula (cf. iDopita et al.l (120001 )). We can identify three 
plausible sources for shock excitation: (i) individual SNe in the stellar cluster, (m) the 
cumulative stellar wind in the cluster, and {iii) the ULX source. Here we make only rough 
estimates of the possible contributions from these three sources; we lack sufficient knowledge 
of the physical conditions in the nebula (e.g., shock velocity, pressure, density, temperature, 
and pre-shock ionization) to apply a full shock-excitation model for the emission lines. 



tc ~ 



The Starburst99 models predict the total supernova production rate in a cluster of age 
4.25 Myr and mass M ~ 10^ M0 and a Salpe ter IMF is about 6 x 10"^ yr ^. Since 



SNRs are visible in emission lines for up to 10^ years (jCioffi et al.lll988l ). the effective number 
of SNRs contributing to the spectrum is at most ~ 6. Typical Ha an d [S ll] luminosities 

[SII] 



from an individual SNR are Ly{n, ~ L\sii] ~ (3 — 6) x 10'^^ ergs s ^ ( iBraun fc Walterbos 



19931 : iMatonick &: FesenI 119971 ). The expected luminosity in these lines from 6 SNRs is then 
~ (1.8 — 3.6) X 10^^ ergs s~^, or 5 to 10 times less than the residual luminosities listed in 
Table 1. 

At the estimated age of 4.25 Myr, the Starburst99 model predicts equal contributions 
from SNe and stellar winds to the mechanical luminosity of the cluster. If the conversion 
efficiencies from mechanical luminosity input to emission line radiation are the same for SNe 
and stellar winds, winds will contribute a comparable amount to the shock-excited lines. 
The combination of SNe and stellar winds still leaves the [S ll] AA6717, 6731 production 
about 2.5 to 5 times below that needed to explain the observed luminosity. 

The third plausible source for shock excitation is the ULX. The presence of such an 
object in the cluster may contribute to the residual line emissions if mechanical energy from 
the X-ray source can be imparted to its surroundin gs in the form o f an outffow. It has 
been postulated that some ULXs are microquasars (IKing et al.l l200ll ). whose radiation is 
coUimated along the jet axis, or supercritical accretion disks with outffows like that in SS433 
(iFabrika fc MeschervakovlboOlh . 
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In the case of microquasars, outflow in the form of a leptonic jet could occur during the 
low state and during outbursts associated with transitions within the very- high/intermediate 
states (Fender, Pooley, & Gallo 2004; Corbel et al. 2003). However, utilizing the known 
correlation between X-ray and radio luminosities for black holes in the low-hard X-ray state 
(Corbel et al. 2003; Gallo, Fender, & Pooley 2003) to estimate the mechanical power in 
the jet and assuming the Chandra observation of CXOU J223706 was obtained when the 
source was near Eddington, then the power in the jet is Lj < 10^^ ergs s~^ even allowing 
for additional transient ejection events (Fender, Belloni, & Gallo 2004). The mechanical 
power imparted to the surrounding nebula in the form of outflow from microquasar-like jets 
from the ULX is then far less than needed to explain the residual line emission. [Note that 
the Galactic microquasars like IE 1740.7-2942 and GRS 1758-258 (Mirabel et al. 1993), for 
instance, have presistent lobe-like emission but are old stellar systems lacking associations 
with young energetic nebulae.] 

In the case of SS433-like objects, outflow is in the form of persisent heavy jets that impart 
a large fraction of their accretion energy into their surroundings in the form of mechanical 
energy. The o bserved kinetic luminosity in the two jets of SS433 is ~ 2 x 10^^ ergs s~^ 



( iFabrikal |200J) . This is comparable to the predicted mechanical luminosity from stellar 
winds and SNe in the star cluster associated with CXOU J223706, 3.5 x 10^^ ergs s~^, at its 
age of 4.25 Myr. 

Likely all three sources of shock excitation contribute to the observed line emission. The 
ULX can be an important, and even dominant, source of radiative and mechanical feedback 
in the overall energetics of the emission region. 



5. Discussion 

We have presented optical photometric and spectroscopic data that help to describe the 
physical environment surrounding the ULX CXOU J223706 in NGC 7331. CXOU J223706 is 
spatially-coincident (within astrometrical errors) with a bright, young, tc = 4.25 ±0.25 Myr, 
star cluster of mass M = (1.1 ± 0.2) x 10^ M0. If the ULX is a part of this bright stellar 
cluster, then the progenitor of the compact accretor must have been at least as massive as 
40-50 Mq in order to have already evolved through the supernova stage. The companion star 
is also likely massive and evolved in order to provide the high mass accretion rate needed to 
power the ULX. 

The emission line spectrum of the surrounding nebula can be interpreted as a result 
of photoinization produced by the cluster UV spectrum, but an additional source of shock 
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excitation is needed to understand the observed [S ll], [Ol] and [Nil] line intensities. The 
mechanical input from SNe and stellar winds from the star cluster can account for only 
~20-40% of the required shock excitation so some additional energy source is needed. The 
obvious candidate is the ULX itself. 

Emission nebulae associated with nearby ULXs are apparently a common phenomenon. 
Although no formal census has been undertaken, we find potential associations with bright, 
extended optical structures for 67 of the 154 ULX candidates tabulated in the ULX survey of 
Swartz et al. (2004). Although these associations are based mainly upon Digitized Sky Survey 
images and must await more conclusive analysis, ma ny detailed studies of ULXs with associ- 



ated nebulae have appeared in the recent li t erature (IPakuU fc Mirioni 



2005: Kuntz et al. 2005; Pakufll et al. 2006; Grise et al 



2006a 



2003 



Lehmann et al. 



bl: lAbolmasov et al.ll2007h 



In particular, the nebula associated with CXOU J223706 in NGC 733 1 can be compare d 
with the compact nebula MF 16 associated with the ULX NGC 6946 X-1 Jsiair et al.lboOlj ). 
The MF 16 nebula is an isolated object: It does not appear to be conta minated by associ- 



ated star clusters or ne arby HII regions. The line luminosities of MF 16 (IBlair et al.l 12001 



Abolmasov et al.ll2006l ) are L^n] ~ ~ 2 x 10^^ ergs s~\ Lh/j ^ Q x 10^^ ergs s"^; equal 
to the residual line luminosities in the NGC 7331 nebula (Table 1). On the other hand, 
[O III] AA4959, 5007 and [O l] A6300 A lines in MF 16 are 3.5 times brighter than those in the 
NGC 7331 residuals. These two lines are excited in different ways: by UV photoionization 
in the case of [O ill] AA4959, 5007 and by shock excitation in the case of [O l] A6300 A. 

Additional sources of mechanical energy in the form of more or less coUimated winds 
or jets is not unprecedented in ULX nebulae. In the nebulae of ULXs NGC 6946 X-1 and 
Holmberg II X-1 (a nother "isolated" ULX nebula), radia l veloc ity gradients of 50—100 km s~^ 
have been found (ILehmann et al.l l2005l : iFabrika et al.l 120061 ) . It was concluded that the 
nebulae have to be powered by the central ULXs and that, therefore, the ULXs must also 
provide a source of gas coUisional excitation in the ULX-associate d nebulae. It was found 
recently through optical spectroscopy of 8 ULX-associated nebulae (lAbolmasov et al.l 120071 ) 
that in at least half of them shock excitation dominates. We cannot say definitely about the 
kinematical properties of the gas in the CXOU J223706 counterpart - we lack both spectral 
and spatial resolution for this - but it is likely the ULX plays an important role here as well. 
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Table 1: Integrated Emission Line Fluxes and Luminosities 
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Fig. 1. — Isaac Newton Telescope Ha image of a portion of the NGC 7331 field showing the 
location of three ULXs in relation to the central regions of the galaxy. At the distance of 
14.6 Mpc to NGC 7331, 1' corresponds to 4.25 kpc. Circles denote 2" radius locations of the 
three ULX candidates detected. CXOU J223706 is the middle source. 
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Fig. 2.— 8"x8" Hubble F450W image of the field around the ULX source CXOU J223706. 
The X-ray source position is marked by the " X" of ±0'.'5 radius corresponding to the combined 
Chandra and Hubble positional uncertainty. Ground-based Ha (black) and [O ill] A5007 
(white) isophotes at the 20, 50 and 85% levels above their backgrounds are superimposed 
(see § 12. 3p . The slit position of the optical spectrum (§ |3]) is also shown. 
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Fig. 3. — Integral spectrum of the cluster near CXOU J223706 (top). The lower thick solid 
curve is the best-fitting stellar spectrum model estimate, divided by a factor of 3 for clarity, 
and the dotted curve is the combination of this model and the best-fitting photoionization 
model (see §11]), divided by a factor of 2. The brightest emission lines including Wolf-Rayet 
blends are superscribed. 



